A meta-analysis of the magnetic line broadening in the solar atmosphere

A multi-line Bayesian analysis of the Zeeman broadening in the solar atmosphere is presented. A hierarchical probabilistic model, based on the simple but realistic Milne-Eddington approximation to the solution of the radiative transfer equation, is used to explain the data in the optical and near infrared. Our method makes use of the full line profiles of a more than 500 spectral lines from 4000 $\AA$ to 1.8 $\mu$m. Although the problem suffers from a strong degeneracy between the magnetic broadening and any other remaining broadening mechanism, the hierarchical model allows to isolate the magnetic contribution with reliability. We obtain the cumulative distribution function for the field strength and use it to put reliable upper limits to the unresolved magnetic field strength in the solar atmosphere. The field is below 160-180 G with 90% probability.Comment: 9 pages, 6 figures, accepted for publication in A&A. Fixed reference

Determination of the cross-field density structuring in coronal waveguides using the damping of transverse waves

Time and spatial damping of transverse magnetohydrodynamic (MHD) kink oscillations is a source of information on the cross-field variation of the plasma density in coronal waveguides. We show that a probabilistic approach to the problem of determining the density structuring from the observed damping of transverse oscillations enables us to obtain information on the two parameters that characterise the cross-field density profile. The inference is performed by computing the marginal posterior distributions for density contrast and transverse inhomo- geneity length-scale using Bayesian analysis and damping ratios for transverse oscillations under the assumption that damping is produced by resonant absorption. The obtained distributions show that, for damping times of a few oscillatory periods, low density contrasts and short inho- mogeneity length scales are more plausible in explaining observations. This means that valuable information on the cross-field density profile can be obtained even if the inversion problem, with two unknowns and one observable, is a mathematically ill-posed problem.Comment: 5 pages, 3 figures, accepte

Bayesian least squares deconvolution

Aims. To develop a fully Bayesian least squares deconvolution (LSD) that can be applied to the reliable detection of magnetic signals in noise-limited stellar spectropolarimetric observations using multiline techniques. Methods. We consider LSD under the Bayesian framework and we introduce a flexible Gaussian Process (GP) prior for the LSD profile. This prior allows the result to automatically adapt to the presence of signal. We exploit several linear algebra identities to accelerate the calculations. The final algorithm can deal with thousands of spectral lines in a few seconds. Results. We demonstrate the reliability of the method with synthetic experiments and we apply it to real spectropolarimetric observations of magnetic stars. We are able to recover the magnetic signals using a small number of spectral lines, together with the uncertainty at each velocity bin. This allows the user to consider if the detected signal is reliable. The code to compute the Bayesian LSD profile is freely available.Comment: 8 pages, accepted for publication in A&

Image Reconstruction with Analytical Point Spread Functions

The image degradation produced by atmospheric turbulence and optical aberrations is usually alleviated using post-facto image reconstruction techniques, even when observing with adaptive optics systems. These techniques rely on the development of the wavefront using Zernike functions and the non-linear optimization of a certain metric. The resulting optimization procedure is computationally heavy. Our aim is to alleviate this computationally burden. To this aim, we generalize the recently developed extended Zernike-Nijboer theory to carry out the analytical integration of the Fresnel integral and present a natural basis set for the development of the point spread function in case the wavefront is described using Zernike functions. We present a linear expansion of the point spread function in terms of analytic functions which, additionally, takes defocusing into account in a natural way. This expansion is used to develop a very fast phase-diversity reconstruction technique which is demonstrated through some applications. This suggest that the linear expansion of the point spread function can be applied to accelerate other reconstruction techniques in use presently and based on blind deconvolution.Comment: 10 pages, 4 figures, accepted for publication in Astronomy & Astrophysic

On the inversion of Stokes profiles with local stray-light contamination

Obtaining the magnetic properties of non-resolved structures in the solar photosphere is always challenging and problems arise because the inversion is carried out through the numerical minimization of a merit function that depends on the proposed model. We investigate the reliability of inversions in which the stray-light contamination is obtained from the same observations as a local average. In this case, we show that it is fundamental to include the covariance between the observed Stokes profiles and the stray-light contamination. The ensuing modified merit function of the inversion process penalizes large stray-light contaminations simply because of the presence of positive correlations between the observables and the stray-light, fundamentally produced by spatially variable systematics. We caution that using the wrong merit function, artificially large stray-light contaminations might be inferred. Since this effect disappears if the stray-light contamination is obtained as an average over the full field-of-view, we recommend to take into account stray-light contamination using a global approach.Comment: 5 pages, 3 figures, accepted for publication in Ap

Compressive Sensing for Spectroscopy and Polarimetry

We demonstrate through numerical simulations with real data the feasibility of using compressive sensing techniques for the acquisition of spectro-polarimetric data. This allows us to combine the measurement and the compression process into one consistent framework. Signals are recovered thanks to a sparse reconstruction scheme from projections of the signal of interest onto appropriately chosen vectors, typically noise-like vectors. The compressibility properties of spectral lines are analyzed in detail. The results shown in this paper demonstrate that, thanks to the compressibility properties of spectral lines, it is feasible to reconstruct the signals using only a small fraction of the information that is measured nowadays. We investigate in depth the quality of the reconstruction as a function of the amount of data measured and the influence of noise. This change of paradigm also allows us to define new instrumental strategies and to propose modifications to existing instruments in order to take advantage of compressive sensing techniques.Comment: 11 pages, 9 figures, accepted for publication in A&